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Computer  Method- 
Obtaining  Radiant  Interchange  Factors  for  a  Quadruped 

D.  G.  Stevens,1  J.  R.  Bohart,!and  LeRoy  Hahn! 


Introduction 


The  exchange  of  an  animal's  heat  to  and  from 
its  surroundings  is  necessary  for  its  life.  Radiation 
is  one  of  the  modes  of  heat  exchange  which  trans- 
fers heat  from  an  animal's  surface  to  its  surround- 
ings. Estimates  of  the  net  radiant  energy  exchange 
for  a  body  require  knowledge  of  the  radiant  inter- 
change factor  (shape  factor).  A  generalized  method 


for  determining  various  radiant  interchange  fac- 
tors for  housed  quadrupeds  (four-legged)  using  a 
computer  program  for  evaluation  is  described  in 
this  report.  The  method  is  applicable  to  any  animal 
housing  situation  in  which  the  surfaces  surround- 
ing the  animal  can  be  assumed  parallel  or  perpen- 
dicular to  the  centerline  of  the  animal. 


Quadruped  and  Enclosure 


A  quadruped  model  is  approximated  by  using 
seven  right  circular  cylinders  (fig.  1).  These  cylin- 
ders represent  the  four  legs,  the  torso,  the  neck 
and  the  head.  The  quadruped  model  assumes  the 
neck  parallel  to  the  torso  and  the  head  at  right 
angles  to  the  neck. 

The  quadruped  is  enclosed  within  six  rectangu- 

Determining  Radiant 

The  radiant  energy  exchange  solution  requires 
finding  the  interchange  factors  between  each  cylin- 
der of  the  model  and  a  particular  enclosing  sur- 
face. The  two  basic  interchange  factors,  a  finite 
cylinder  to  a  triangle  and  a  circle  to  a  triangle, 
were  obtained  from  Tripp,  Hwang  and  Crank.3  The 
position  of  the  triangle  in  relation  to  the  cylinder 
and  circle,  the  equations  and  the  subroutines  to 
solve  these  equations  are  found  in  the  appendix. 
These  integral  equations  were  solved  numerically 


'Agricultural  engineers,  bioengineering  research  unit, 
Agricultural  Research  Service,  U.S.  Department  of  Agri- 
culture, Columbia,  Mo.  65201 

2 Mechanical  Engineer,  Babcock  &  Wilcox  Co.r  Power 
Generation  Div.,  Lynchburg,  Va.  24505 

3Tripp,  Wilson,  Ching-Lai  Hwang,  and  R.  E.  Crank. 
Radiation  Shape  Factors  for  Plane  Surfaces  and  Spheres, 
Circles  or  Cylinders.  Kans.  State  Univ.  Bui.  Spec.  Rpt. 
No.  16,  vol.  46,  No.  4.    April  1962. 


lar  surfaces,  either  parallel  or  perpendicular  to  the 
centerline  of  the  cylinders.  A  generalized  computer 
program  was  developed  using  as  input  variables 
the  size  of  the  walls  and  the  quadruped's  com- 
ponent cylinder  and  position  of  the  animal  within 
the  enclosure.  This  allows  use  of  the  program  with 
many  enclosure  configurations  and  various  sizes 
of  quadruped. 

Interchange  Factors 

and  generalized  in  a  computer  program.  The  re- 
maining shape  factors  around  the  sides  and  ends 
of  the  component  cylinders  were  obtained  by  flux 
algebra,  using  methods  described  by  Wiebelt," 
in  the  main  computer  program  or  its  subroutines.5 
The  list  of  assumptions  made  in  writing  the 
program  follows: 

1.  The  legs  cannot  radiate  to  the  ceiling. 

2.  The  head  and  neck  do  not  radiate  to  the 
rear  walls. 

3.  Only  the  lower  section  of  the  head  (below 
the  lower  side  of  the  neck  cylinder)  radiates  to 
the  sidewalls. 


"Wiebelt,  J.  A.  Engineering  Radiation  Heat  Transfer. 
Holt,  Rinehart  and  Winston,  New  York.  1966. 

5A  description  and  program  listing  of  the  main  program 
and  subroutines  (FORTRAN)  as  developed  for  use  on  the 
IBM  370/165  computer  may  be  obtained  by  writing  the 
Missouri  Climatic  Laboratory,  Agr.  Engin.  Dept.,  Univ.  of 
Mo.,  Columbia,  Mo.  65201 
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Figure  1.— Placement  of  cylinders  to  form  model 
quadruped. 


4.  The  neck  cylinder  radiates  to  the  section  of 
the  radiation  panels  that  extends  in  front  of  the 
body  cylinder. 

5.  The  quadruped  is  standing  centered  between 
the  sidewalls  of  the  chamber. 

6.  The  required  shape  factors  are  calculated  for 
an  application  only  in  which  all  surfaces  have  black 
body  characteristics. 

Example  Computations 

As  an  example  application  to  a  specific  case,  the 
program  was  used  for  a  quadruped  representing  a 
Holstein  cow  with  dimensions  approximating  those 
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given  by  Brody.6  The  radius  and  length  (centi- 
meters) of  the  cylinders  used  are: 

Radius  Length 

Legs    3.15  89.0 

Body  32.00  127.5 

Neck  15.25  57.4 

Head  10.16  54.6 

The  enclosure  in  this  case  is  an  individual  envi- 
ronmental chamber  in  the  Missouri  Climatic  Lab- 
oratory. The  inside  dimensions  of  the  chamber  are 
285  cm  long  x  122  cm  wide  x  193  cm  high.  The 
sidewalls  have  three  radiation  panels  and  the  ceil- 
ing has  two  radiation  panels  placed  as  shown  in 
figures  2  and  3. 

By  definition,  the  sum  of  the  shape  factors  from 
one  surface  to  all  surrounding  surfaces  must  equal 
1.0.  The  shape  factors  for  the  body  cylinder  and 
rear  body  circle  were  used  to  check  the  validity  of 
the  program.  The  error,  represented  by  the  devia- 
tion from  unity,  for  the  sum  of  the  body  cylinder 
shape  factors  was  0.132  percent  and  for  the  rear 
body  circle  it  was  0.44  percent. 

Shown  below  is  an  example  for  finding  the  shape 
factor  from  the  left  rear  leg  to  the  rear  wall.  Five 
steps  are  used: 

Step  1.  The  shape  factor  from  the  legs  to  that 
portion  of  the  wall  that  is  below  the  top  of  the  leg 
cylinder  and  to  the  right  of  the  leg  centerline.  This 
uses  the  subroutine  FAC3.  (Shown  in  Case  II  in 
the  appendix.) 

Step  2:  This  step  is  the  same  as  Step  1  but 

6 Brody,  Samuel.  Bioenergetics  and  Growth.  Reinhold 
Publishing  Co.,  New  York,  p.  650.  1945. 
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Figure  2. — Ceiling  of  individual  cow  chambers. 


Figure  3. — Wall  of  individual  cow  chambers. 


finds  the  shape  factor  from  the  leg  to  portion  of 
the  wall  to  the  left  of  the  leg  centerline. 

Step  3:  The  shape  factor  from  the  leg  to  that 
portion  of  the  wall  that  is  above  the  top  of  the  leg 
cylinder  and  to  the  right  of  the  leg  centerline.  This 
uses  the  subroutine  FAC8  (shown  in  Case  I  in  the 
appendix). 

Step  4:  This  step  is  the  same  as  Step  3  but  to 
the  portion  of  the  wall  to  the  left  of  the  leg 
centerline. 

Step  5:  The  four  shape  factors  found  in  the 


preceding  steps  are  summed  to  find  the  shape 
factor  from  the  left  leg  to  the  rear  wall. 

The  remaining  shape  factors  are  found  by  apply- 
ing the  two  subroutines  in  various  combinations. 
The  complete  set  of  shape  factors  found  by  the 
main  program  is  shown  in  table  1. 

Refinements  and  Application 

The  initial  program  computes  shape  factors  that 
must  be  considered  a  first  approximation  because 
of  the  number  of  assumptions  that  were  made. 


Table  1. — Radiation  shape  factors  for  all  components 


Emitting  surface 
Right  rear  leg  


Right  front  leg 


Left  rear  leg 


Left  front  leg 


Body  (torso) 


Body  receiving  surface  Shape  factor 

Right  lower  radiation  panel  0.136387 

Right  middle  radiation  panel  .059776 

Right  upper  radiation  panel  .007587 

Left  lower  radiation  panel  .041207 

Right  wall  exclusive  of  radiation  panels  .091861 

Rear  wall  .150899 

Floor  .277738 

.  Right  lower  radiation  panel  .158680 

Right  middle  radiation  panel  .068826 

Right  upper  radiation  panel  .008038 

Left  lower  radiation  panel  .049300 

Right  wall  exclusive  of  radiation  panels  .062422 

Front  wall  .143710 

Floor  .278445 

.  Right  lower  radiation  panel  .041207 

Left  lower  radiation  panel  .136387 

Left  middle  radiation  panel  .059776 

Left  upper  radiation  panel  .007587 

Left  wall  exclusive  of  radiation  panels  .091861 

Rear  wall  .150899 

Floor  .277738 

.  Right  lower  radiation  panel  .049300 

Left  lower  radiation  panel  .158680 

Left  middle  radiation  panel  .068826 

Left  upper  radiation  panel  .008038 

Left  wall  exclusive  of  radiation  panels  .062422 

Front  wall  .143710 

Floor  .278445 

.  Right  lower  radiation  panel  .052813 

Right  middle  radiation  panel  .091062 

Right  upper  radiation  panel  .075321 

Right  ceiling  radiation  panel  .078737 

Left  lower  radiation  panel  .052813 

Left  middle  radiation  panel  .091062 

Left  upper  radiation  panel  .075321 

Left  ceiling  radiation  panel  .078737 

Right  wall  exclusive  of  radiation  panels  .086526 

Left  wall  exclusive  of  radiation  panels  .086526 

Ceiling  exclusive  of  radiation  panels  .041240 

Front  wall  .016158 

Rear  wall  .017494 

Floor  .154872 

Total  (body)  0.998678 
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Table  1. — Radiation  shape  factors  for  all  components — Continued 


Emitting  surface 


Rear  body  circle  receiving  surface 


Shape  factor 


Neck   Left  or  right  lower  radiation  panel  0.003462 

Left  or  right  middle  radiation  panel  .032240 

Left  or  right  upper  radiation  panel  .028295 

Left  or  right  ceiling  radiation  panel  .063411 

Left  or  right  wall  exclusive  of  radiation  panels  .166015 

Ceiling  exclusive  of  radiation  panels  .122637 

Head   Left  or  right  lower  radiation  panel  .009473 

Left  or  right  middle  radiation  panel  .054247 

Left  or  right  upper  radiation  panel  .010922 

Left  or  right  wall  exclusive  of  radiation  panels  .062841 

Front  wall  .412551 

Front  body  circle  ....  Left  or  right  lower  radiation  panel  .007960 

Left  or  right  middle  radiation  panel  .030836 

Left  or  right  upper  radiation  panel  .016153 

Left  or  right  ceiling  radiation  panel  .009692 

Rear  body  circle  Right  lower  radiation  panel  .001747 

Right  middle  radiation  panel  .008116 

Right  upper  radiation  panel  .003331 

Right  ceiling  radiation  panel  .002854 

Left  lower  radiation  panel  .001747 

Left  middle  radiation  panel  .008116 

Left  upper  radiation  panel  .003331 

Left  ceiling  radiation  panel  .002854 

Right  wall  exclusive  of  radiation  panels  .143718 

Left  wall  exclusive  of  radiation  panels  .143718 

Ceiling  exclusive  of  radiation  panels  .085124 

Rear  wall  .550211 

Floor  .050806 


Total  (rear  body  circle) 


1.005672 


Refinements  that  could  be  incorporated  are  ap- 
proximating the  quadruped  using  right  circular 
ellipses,  truncated  cylindrical  cones  or  truncated 
elliptical  cones. 

Physical  models  of  total  heat  balance  for  quad- 


rupeds as  developed  for  the  pig  by  Beckett  and 
Vidrine7  require  evaluation  of  radiation  heat  trans- 
fer. The  method  of  obtaining  radiant  interchange 
factors  outlined  in  this  report  should  expedite 
estimating  the  radiation  component  of  the  total 
heat  balance. 


The  shape  factor  for  radiation  from  a  finite  cyl- 
inder to  a  surface  that  is  parallel  to  the  centerline 
of  the  cylinder  cannot  be  determined  directly.  The 
radiation  being  received  by  surface  B,  in  figure  4, 
and  the  two  right  triangles  connecting  the  surface 
and  the  cylinder  is  equal  to  the  radiation  leaving 
the  shaded  portion,  A,  of  the  cylinder.  Therefore, 
to  find  the  radiation  leaving  the  cylinder  and  strik- 
ing the  surface  B,  that  is,  the  shape  factor,  the 
radiation  passing  through  the  two  triangles  is  sub- 
tracted from  that  leaving  A. 


Finite  Cylinder  To  a  Right  Triangle 

Figure  5  shows  a  view  of  the  cylinder  and  the 
placement  of  the  right  triangle  for  development  of 
the  equation.  The  detailed  explanation  (next  page) 
of  the  solution  of  this  shape  factor  is  given  by 
Tripp  and  others.8 


'Beckett,  F.  E.  and  Clyde  Vidrine.  A  Mathematical 
Model  of  Heat  Transfer  in  a  Pig.  Trans,  of  ASAE  13(4): 
450-454.  1970. 

'  See  footnote  3,  page  1. 
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Figure  4. — Finite  cylinder  to  a  parallel  surface. 


Figure  5. — Placement  of  right  triangle  to  develop  equation. 


2  1 

p  =  p      s       p  = 

Al>A2        12      Al    A2*A1  2-n2W 


1 


V777 


tan 


+  (q+1) (q-1) 


VU2+(a+l)2}  |f^+(a-l)2} 

-/ 


c2 


tan1  J 


(a-1)  {  ft^+to+l)2} 
(a+1)  W2+(q-l)2} 


ada 


V2  +  (q+1) (q-1) 


/(q-1)  {ft^+fa+D2} 

V{f/2+(q+l)2}     |^2+(q-l)2}  f  {(^+(q-l)2} 


tan 


-la  , 
sec       —  ada 

c 


where 


B  =  —t  C  =  —,  W  =  — ,  q  =  — ,  Is  an  increment  of 

K  n  n  n 


This  equation  is  evaluated  in  the  subroutine 
FACTOR  (next  page)  which  uses  numerical  quad- 
rature to  evaluate  the  integrals. 
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SUBROUTINE  FACTOR  (ALl,AL2,H,R, AAREA) 


C  C 

C***** ****************************************************** ****************£ 

c  c 

C      THIS  SUBROUTINE  FINDS  THE  SHAPE  FACTOR  FROM  A  CYLINDER  TO  A  RIGHT  C 

C      TRIANGLE  WHICH  IS  PERPENDICULAR  TO  THE  CENTERLINE  AND  EVEN  WITH  ONE  C 

C      END  OF  THE  CYLINDER.  C 

£              AL1      =  LENGTH  OF  OPPOSITE  SIDE  OF  TRIANGLE  £ 

£              AL2          LENGTH  OF  ADJACENT  SIDE  OF  TRIANGLE  Z, 

H              HEIGHT  OF  CYLINDER  £ 

Z,  R  RADIUS  OF  CYLINDER  £ 
p              AAREA  =  SHAPE  FACTOR 

Q* *************************************************************** ***********C 

C  C 


DOUBLE  PRECISION  A , ZA ,B , C ,W,W2 ,AP1 ,AM1 ,TBC ,AR(309) ,B2 ,C2 

DOUBLE  PRECISION  AP12 ,AM12 , ZB1 , ZB2 , ZB ,ZC , ZD ,ANC 

DOUBLE  PRECISION  AREA1 ,AREA2 

B=AL1/R 

C=AL2/R 

W=H/R 

ITTER=300 

B2=B**2 

C2=C**2 

ANC= (DSQRT (B2+C2) -1 . 0) /FLOAT (ITTER) 

IT=ITTER+1 

A=l. 000001 

W2=W**2 

TBC=DATAN(B/C) 

DO  1  1=1, IT 

AP1=A+1 . 0 

AM1=A-1.0 

ZA=DATAN (DSQRT (API /AMI) ) 

AP12=AP1**2 

AM12=AMl**2 

ZB1=W2+AP12 

ZB2=W2-hAM12 

ZB=W24AP1*AM1 

ZC=AB/DSQRT(ZB1*ZB2) 

ZD=DATAN (DSQRT ( (AM1*ZB1) / (AP1*ZB2) ) ) 
AR(I)=TBC*(ZA-ZC*ZD)*A 
1  A=A+ANC 

CALL  QUADR (AR, IT ,ANC, AREA 1, ITER) 
ANC= (DSQRT (B2+C2) -C) /FLOAT (ITTER) 
A=C 

DO  3  1=1, IT 

AP1=A+1 . 0 

AM1=A-1.0 

AP12=AP1**2 

AM12=AM1**2 

ZB1=W2+AP12 

ZB2=W24AM12 

ZB=W2+AP1*AM1 

ZC=ZB/DSQRT (ZB1*ZB2) 

ZD=DATAN (DSQRT ( (AM1*ZB1) / (APl*ZB2) ) ) 
AR(I)=DARCOS (C/A)*A*(ZA-ZC*ZD) 
3  A=A+ANC 

CALL  QUADR  (AR,  IT  ,ANC,AR£A2,  ITER) 
AAREA= (0 . 0506606/W) *(AREA1-AREA2) 
RETURN 
END 
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c  c 

c  c 

CSUBROUTINE      QUADRATURE  C 

C  C 

C  PURPOSE:  C 
C          INTEGRATES  A  GIVEN  TABULATED  FUNCTION  AT  A  SET  OF  EQUALLY  SPACED  C 

C          POINTS .  C 

C  C 

C    USAGE :  C 

C          CALL  QUADR(Z,N,H,S  IER)  C 

C  C 

C    DESCRIPTION  OF  PARAMETERS :  C 

C          Z      =  A  VECTOR  OF  LENGTH  N  CONTAINING  THE  VALUES  OF  THE  FUNCTION  C 

C                     TO  BE  INTEGRATED  C 

C          N          THE  NUMBER  OF  FUNCTION  VALUES  TO  BE  INTEGRATED  C 

C          H          THE  SPACING  OF  THE  FUNCTION  VALUES  C 

C  S  THE  RESULTANT  VALUE  OF  THE  COMPUTED  INTEGRAL  C 
C          IER  =  RESULTANT  ERROR  CODE  WHERE  IER=0  INDICATES  NO  ERROR,  IER=1  C 

C                     INDICATES  N  LESS  THAN  2,  IER=2  INDICATES  H=0 .  C 

C  C 

C    REMARKS:  C 

C          SEE  ERROR  CODES  ABOVE.  C 

C  C 

C    METHOD :  C 

C          NUMERICAL  QUADRATURE  IS  PERFORMED  USING  A  5  POINT  FORMULA.  UNTIL  C 

C          FEWER  THAN  5  POINTS  REMAIN.    A  4,  5,  OR  2  POINT  FORMULA  IS  THEN  C 

C         USED  TO  COMPLETE  THE  INTEGRATION.  C 

C  C 

Q  'k^'k''k~k~k~k~k  "k~k  tWt  "k^k  itirk'k'k  irfrk  k  "k"k"k-Jck  'Jcfrfckicfc  "k~fc^dcfck'7C^ck~k  'k'k'kick  irk  irkifick^ckickrk'kirk'k'k'k'k  ~k  "kQ 

c  c 

SUBROUTINE  QUADR(Z,N,H,S  ,IER) 

c  c 

c  c 

C  IF  A  DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED,  THE  C  IN  C 
C    COLUMN  1  SHOULD  BE  REMOVED  FROM  THE  DOUBLE  PRECISION  STATEMENT  WHICHC 

C    FOLLOWS.  C 

c  c 

Q  ick'k'irkk'k'k'k  fddCft  "k"k~k~k~k"k  ~k~k~k~k^<~k"k      ~k~k~k  ~^rk~k  tWt  ~k~k~k~k  ~k~k~k~k  "kirk~k  ~k~k~k~k~kk  "k~k~k~k~k~k~k  ~k~k  mk~k~k~k~k~k  A£ 

c  c 

DOUBLE  PRECISION  Z,H,S 

s=.o 

IF (N-l)32, 32,33 

32  IER=1 
RETURN 

33  IF(H)34,35,34 
35  IER=2 

RETURN 

34  IF(N-5)40,38,38 

C  C 

C         CALCULATE  USING  5  POINT  FORMULA  C 

C  C 

38  DO  39  1=5, N,4 

39  S=S+7.*Z(I-4)+32 .*Z(I-3)+12 .*Z(I-2)+32 .*Z(I-l)+7  .*Z(I) 
S=S*2 . /45 . 

40  J=N-(N/4)*4+l 

GO  TO  (45, 50, 47, 48), J 

C  C 

C        CALCULATE  USING  4  POINT  FORMULA  C 

C  C 
45  S=S+.375*(Z(N-3)+3.*Z(N-2)+3.*Z(N-l)+Z(N)) 

GO  TO  50 

C  C 

C        CALCULATE  USING  2  POINT  FORMULA  C 

c  c 

47  S=S+(Z(N-l)+Z(N))/2. 
GO  TO  50 

c  c 

C        CALCULATE  USING  3  POINT  FORMULA  C 

C  C 

48  S=S+(Z(N-2)+4.*Z(N-l)+Z(N))/3. 
50  S=S*H 

IER=0 

RETURN 

END 
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Finite  Circle  To  a  Right  Triangle 

The  shape  factor  found  from  this  subroutine  is 
used  to  find  the  shape  factors  from  the  end  circles 
of  the  body  to  the  walls,  shown  in  figure  6.  A 
detailed  explanation  of  the  procedure  is  given  by 
Tripp  and  others' 

-1  z 
tan  J 

-  j  (Z2  +  rV  +  sec2*)2  -  4X222  sec2(M*  j 


The  preceding  equation  is  evaluated  in  the  sub- 
routine CIRFAC. 


See  footnote  3,  page  1. 


Figure  6 — Placement  of  the  right  triangle  in  relation  to 
the  circle  and  the  numbers  used  in  the  equation  to  find 
the  shape  factor. 


SUBROUTINE  CIRFAC  (ALl  ,AL2  ,D,R, AAREA) 

C 

c 

THIS  SUBROUTINE  FINDS  THE  RADIATION  SHAPE  FACTOR  FROM  A  CIRCLE  TO  A 
RIGHT  TRIANGLE  THAT  IS  PARALLEL  TO  THE  CIRCLE. 


C 
C 

c 
c 
c 
c 
c 
c 
c 

c 


AL2 
ALl 
D 
R 


=  LENGTH  OF  ADJACENT  SIDE  OF  TRIANGLE 

=  LENGTH  OF  OPPOSITE  SIDE  OF  TRIANGLE 

=  NORMAL  DISTANCE  BETWEEN  CIRCLE  AND  TRIANGLE 

=  RADIUS  OF  CIRCLE 


AAREA  =  SHAPE  FACTOR 


C 
C 

c 

c 
c 
c. 
c 
c 
c 
c 

fciricfcicft  fcfrtcicft  A"A-  'A™ fcfciriricyririric^ricfc  A'A-A-^A  ".Wr  sfcfe  ^ficfcSddclFx  xxx  A-  A*  *A"A  A  x  A  A-  ft  Q 

c 

DOUBLE  PRECISION  X,Y  ,Z  ,X2  ,TZY  ,  Z2  ,ANC  ,A , CONST  ,SEC2  ,AR(309)  ,AREA1 

X  =  R/D 

Y  =  D/AL1 

Z  =  D/AL2 

X2  =  X*X 

TZY  =  DATAN  (Z/Y) 

Z2  =  Z*Z 

ITTER  =300 

ANC  =  TZY /FLOAT  (ITTER) 
IT  =  ITTER+1 
A  =  0.0 

CONST  =  (1.0/(Y*Z)+(1.0+X2)*TZY) 
DO  1  1=1, IT 

SEC 2  =  (1.0/DCOS(A))**2.0 

AR(I)  =  (DSQRT((Z2+X2*Z2+SEC2)**2.0-4.0*X2*Z2*SEC2))/Z2 
1  A  =  A+ANC 

CALL  QUADR  (AR ,  IT  ,ANC  ,AREA1 ,  ITER) 
AAREA  =  (CONST-AREA1)/(4.0*3.1417*X2) 
RETURN 
END 
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Circle  To  a  Rectangular  Wall 

In  this  subroutine,  the  wall  is  perpendicular  to 
the  plane  of  the  circle,  shown  in  figure  7. 
The  shape  factor  ^i-*"1 
is  found  from  the  following  equations: 


■i  "  Fc 


(1  +  2)  -  'C. 


1      v  "1  , 

CALL  CIRFAC  (A.B,D,R,F1) 


(1  +  2) 


1  -  -1 
-r  tan 


The  CIR1  subroutine  implements  the  preceding 
equations. 


Figure  7. — End  circle  in  relation  to  a  wall. 


SUBROUTINE  CIR1 (A, B ,D , R, F) 


c  c 

c  c 

C  THIS  SUBROUTINE  FINDS  THE  SHAPE  FACTOR  FROM  A  CIRCLE  TO  A  RECTANGLE  C 

C  EVEN  AND  NORMAL  TO  THE  PLANE  OF  THE  CIRCLE.  C 

C  C 

C  A  =  HEIGHT  OF  RECTANGLE  C 

C  B  =  RADIAL  DISTANCE  FROM  CENTER  OF  CIRCLE  TO  RECTANGLE  C 

C  D  =  LENGTH  OF  RECTANGLE  C 

C  R  =  RADIUS  OF  CIRCLE  C 

C  F  =  SHAPE  FACTOR  C 

c  c 

c****** ****************************** ********************************** *q 

c  c 


CALL  CIRFAC (A, B,D,R,F1) 

F2  =  ATAN(A/B)/(2. 0*3. 1417) 

F  =  F2-F1 

RETURN 

END 


Case  I 

In  this  case,  the  radiation  shape  factor  to  be 
found  is  from  a  cylinder  to  a  wall  that  is  parallel  to 
the  centerline  of  the  cylinder  and  the  bottom  of 
the  wall  is  the  top  of  the  cylinder,  from  cylinder  CI 
to  Wall  1. 

The  equation  is  numbered  as  in  figure  8,  and 
the  radiating  shape  factor  Fc1->-i  is  found  from 
the  equations  that  follow. 


P„        ,  =  CALL  FACTOR  (B,A.H,R,F1T4) 

1  A  A 

F  =    C    +  %  Ar 

C  1       2  Ar  C2^3 

''I 

F(C    +  c  )->-3  =  CALL  FACT0R  (B,A,H  +  C,R,F12T3) 
Fc  =  CALL  FACTOR  (B,A.C,R,F2T3) 


Figure  8. — Placement  of  cylinder  CI  and  Wall  1. 
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The  FAC8  subroutine  implements  the  above  formula. 


SUBROUTINE  FAC8(A,B,C,H,R,F) 


c  c 

£*************************************************************************c 

C  THIS  SUBROUTINE  FINDS  THE  RADIATION  SHAPE  FACTOR  FOR  A  CYLINDER  C 

C  TO  A  WALL  WHICH  IS  PARALLEL  TO  THE  CENTERLINE  OF  THE  CYLINDER  AND  C 

C  THE  BOTTOM  OF  THE  WALL  IS  AT  THE  TOP  OF  THE  CYLINDER.  C 

C  C 

C  THE  PARAMETERS  ARE-  C 

C  C 

C              A  =  THE  PERPENDICULAR  DISTANCE  FROM  THE  WALL  TO  THE  CYLINDER  C 

C              B  =  THE  LENGTH  OF  THE  WALL  C 

C              C  =  THE  HEIGHT  OF  THE  WALL  C 

C              H  =  THE  HEIGHT  OF  THE  CYLINDER  C 

C              R  =  THE  RADIUS  OF  THE  CYLINDER  C 

C              F  =  THE  SHAPE  FACTOR  FROM  THE  CYLINDER  TO  THE  WALL  C 

c*************************************************************************c 

c  c 


ANGLE=ATAN(B/A) 
AC 1=R*H*ANGLE 
AC2=R*C*ANGLE 
CALL  FACTOR  (B,A,C,R,F2T3) 
CALL  FACTOR(B,A,H+C,R,F12T3) 
F1T3=( (AC1+AC2) /AC1) *F12T3- (AC2/AC1)*F2T3 
CALL  FACTOR  (B ,A,H, R.F1T4) 
F=F1T4-F1T3 
RETURN 
END 

Case  II 

In  this  case,  the  radiation  shape  factor  to  be 
computed  is  from  a  cylinder  to  a  wall  in  which  the 
wall  is  the  same  height  as  the  cylinder,  and  the 
bottom  of  the  wall  is  even  with  the  bottom  of  the 
cylinder.  The  numbers  assigned  to  the  surfaces  in 
figure  9  are  used  in  the  equations  for  determining 

The  following  are  equations  to  find  F"1~^r 

CALL  FACTOR  (B,A,H,R,F1T1) 
FAC3  subroutine  implements  the  above  equations: 

FC       2  =  FC        (1+2+3)  ~  FCx-*-±  ~  FCx-*~  3 

\+  (1+2+3)  =  h  t3n_1  (BM) 

F„  =  F         -  -  CALL  FACTOR  (B ,A,H ,B,F1T1) 

Ll      1        Cl  J 


Figure  9. — Placement  of  cylinder  Cl  and  Wall  2. 


SUBROUTINE  FAC3 (A,B ,H,R,F) 


c  c 

C* **************************************** A*** ******** ****** ****************c 

C  THIS  SUBROUTINE  FINDS  THE  RADIATION  SHAPE  FACTOR  FROM  A  CYLINDER  C 

C  TO  A  WALL  WHICH  IS  PERPENDICULAR  TO  THE  CYLINDER  AND  THE  SAME  C 

C  HEIGHT .  C 

C  C 

C  THE  PARAMETERS  ARE-  C 

C                 A  =  PERPENDICULAR  DISTANCE  FROM  THE  CYLINDER  TO  THE  WALL  C 

C                 B  -  LENGTH  OF  THE  WALL  C 

C                 H  -    HEIGHT  OF  THE  WALL  AND  CYLINDER  C 

C                 R  =  RADIUS  OF  THE  CYLINDER  C 

C                  F  =  RADIATION  SHAPE  FACTOR  FROM  CYLINDER  TO  WALL  C 

C  C 
c *************** *************************************** *********************c 

c  c 


CALL  FACT0R(B,A,H,R,F1T1) 
AND  -  0.159155  *  ATAN  (B/A) 
F  =  ANG-.2*F1T1 
RETURN 

END  10 


